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An experimental proposal to test the existence of a liquid-liquid critical point in water, based on 
high pressure calorimetric measurements ,is presented on this paper. Considering the existence of 
an intramolecular correlation in the water molecule we show how the response of the specific heat 
at high pressure is different depending on the existence, or not, of the second critical point. If the 
liquid-liquid critical point hypothesis is true there must be a maximum in the specific heat at some 
temperature T > Th for any pressure P > P c (being Th the homogeneous nucleation temperature 
and P c the pressure of the second critical point). This maximum does not appear for the singularity 
free scenario. 



Liquid water exists in a metastable supercooled state 
far bellow the melting temperature. A number of thermo- 
dynamic response functions of supercooled water such as 
the isothermal compressibility and the constant pressure 
specific heat, show a power-law divergence at T=228K 
and P=latm [3,0]. In order to explain these anomalies 
several theories have been proposed. 

The liquid-liquid phase transition hypothesis Q pro- 
poses the existence of a first order line of phase transi- 
tions separating two liquid states of different densities: 
the high density liquid (HDL) and the low density liquid 
(LDL) state. This line of phase transitions has a nega- 
tive slope in the P — T phase diagram and ends up in a 
critical point at T c ~ 200K and P c ~ 1.7kbars. Another 
possible scenario is singularity-free 0]. In this proposal 
the anomalies found in the experiments are not consid- 
ered to end up in real singularities but in non divergent 
maxima due to the anticorrelated fluctuations of volume 
and entropy. 

The direct observation of a possible first order liquid- 
liquid phase transition line and a second critical point in 
water has been hampered by the intervention of the ho- 
mogeneous nucleation process which takes place at higher 
temperatures Tu(P) than the ones corresponding to the 
hypothesized coexistence line. In this paper we are go- 
ing to focus our attention in the following question: Can 
we expect any difference (at temperatures T > Th{P)) in 
the calorimetric response of supercooled water depending 
on the existence, or not, of a second critical point? 

In order to give a possible answer to this question we 
are going to consider a liquid model where the existence, 
or not, of a second liquid-liquid critical point is easily 
tunable. The model is described in detail in [5]. In the 
liquid phase, each water molecule interacts with another 
four water molecules forming hydrogen bonds with en- 
ergy — J. The formation of a hydrogen bond increases 
the volume of the system by a certain amount SV. In 
the simplest version of this model, there is no correla- 
tion between the possible four hydrogen bonds formed 
by a particular water molecule. The correlation length 



related to the formation of the hydrogen bond network 
does not diverge and the thermodynamic response func- 
tions present a simple (non-divergent) maximum at a 
temperature Th(P). This temperature, Th(P), is always 
proportional to the absolute value of the interaction en- 
ergy due to the formation of a hydrogen bond at a given 
pressure, T h (P) oc (J — (PSV)) [|. Note that T h (P) cx J 
for P = and T h (P) = for P = J/SV. The behavior of 
this simple model is the one described by the singularity 
free scenario. 

However by considering an internal correlation among 
the four hydrogen bonds formed by each single molecule 
the scenario changes from singularity- free to critical 
point like This internal degree of correlation in the 
water molecules is tuned by a new intramolecular term 
with energy — J a added to the Hamiltonian of the liquid 
5]. If J a — > oo the water molecule is always internally 
correlated and, for any value of the pressure, the ther- 
modynamic functions are always divergent at Th(P). 

The liquid-liquid critical point at T c and P c shows up 
when J fj is finite and smaller than J. In this case, the wa- 
ter molecules are able to correlate themselves internally 
only if T < T* , being T* a pressure independent tempera- 
ture proportional to J a . For P < P c the water molecules 
are not internally correlated at Th(P) (i.e. Th(P) > T*), 
the correlation length of the hydrogen bonds network is 
unable to tend to infinity, and the response of the ther- 
modynamic functions is non-divergent. 

However, as the pressure increases Ifc(P) decreases. 
So, there is a critical value of the pressure P = P c where 
the water molecules are internally correlated at Th(P c ) 
(i.e. Th(P c ) < T„), the correlation length of the hydro- 
gens bond network is able to go to infinity, and the re- 
sponse of the thermodynamic functions is divergent at 
Th(P c ) = T c . Of course, for any P > P c , the four hy- 
drogen bonds formed by any water molecule are always 
internally ordered at Th (P) < T* and the response of the 
thermodynamic functions are going to be the ones cor- 
responding to a first order phase transition between two 
liquids with different density. 
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FIG. 1: Three possible scenarios depending on the value of 
J CT . (a) If J a = there is no internal correlation in the water 
molecule for any pressure (T* = 0) and there is no long range 
correlation in the formation of the hydrogen bonded network. 
The transition is always singularity- free (dotted line), (b) 
If ,l a < J and different from zero, there is long range cor- 
relation in the transition to a low density liquid only when 
the transition temperature Th is smaller than the tempera- 
ture T* corresponding to the internal correlation of the water 
molecules. In this particular case we obtain a critical point 
and a first order phase transition for any pressure P > P c and 
for any temperature T < T c ~ T* (gray region and continuous 
line), (c) If J a is large enough compared to J there is always 
internal correlation in the water molecules for any pressure 
and the transition to a low density liquid is first order for any 
pressure (gray region and continuous line). 



This mechanism is explained in detail on Fig.l. Note 
how, by tunning the intra-molecular term J a , we can 
easily change the physical mechanism of the simulated 
liquid form singularity free (J CT = 0) to liquid-liquid phase 
transition (0 ^ J a < J). 

In real water, this intra-molecular term in the water 
molecule could be related to the value of the H-O-H 
angle. Experiments show that the relative orientations 
of the arms in the water molecule are correlated, with 
the average H-O-H angle equal to 104.45° in the isolated 
molecule, 104.474° in the gas and 106° in the high T liq- 
uid Q , suggesting an intramolecular interaction between 
the arms. This interaction must be finite because the an- 
gle changes with T, consistent with ab-initio calculations 
Q and molecular dynamics simulations 



Based on the model, we propose the following hypoth- 
esis: If there is a second critical point at P c and we cool 
down our liquid deep into the supercooled region at con- 
stant P > P c , there is going to be internal correlation 
inside the water molecule at a temperature T* > Th (i.e. 
before homogeneous nucleation takes place). The fluc- 
tuations on the energy of the system due to this intra- 
molecular ordering give rise to a maximum in the specific 
heat with the following characteristics: (i) It is located at 
a temperature larger than the one corresponding to nu- 
cleation (T* > Th), (ii) it is non-divergent (because it is 
not due to the cooperative effect of many molecules but 
to the internal ordering of each molecule separately) and 
(iii) the position of the maximum is going to be pressure 
independent (i.e. since J a is constant and does not de- 
pend on pressure the value of T* should be also pressure 
independent). 

On the other hand if P < P c and we cool down our 
liquid, nucleation takes place before intra-molecular or- 
dering (T* < Th) and no maximum appears in the spe- 
cific heat (only the regular nucleation maximum reach- 
able in numerical simulations but unreachable in a real 
supercooled bulk water experiments). In the case of a 
singularity free scenario where there is no critical point 
(critical pressure) , there is never any internal ordering in 
the water molecules and, independently of the pressure 
applied, no maximum in the specif heat is found. 

In order to test this hypothesis we have simulated the 
liquid phase of the model described in Ref. [5J] by com- 
puting the total energy density of states 0. We have 
chosen the following values for the parameters: J/e = 50 
(being e the value of the van der Waals interaction en- 
ergy which in our case, since we are only considering 
the condensed liquid phase, does not play any role) and 
SV/v = 50 (being v a the hard core volume of the wa- 
ter molecule). We have performed two sets of calcula- 
tions, one for a liquid with no liquid-liquid critical point 
(Jff/e = 0) and another one with a liquid-liquid critical 
point (Jtr/e = 5). All other parameters are the same 
than in Ref. [5j. The behavior of energy fluctuations 
(constant pressure specific heat) has been calculated for 
a system of 225 water molecules with periodic boundary 
conditions and for the following values of the pressure 
Pv /e = 0.54,0.74,0.84,0.89,0.94,0.955,0.97. Results 
are shown in Fig. 2 

Fig. 2 (a) are the results for the liquid with no-critical 
point. In this case no secondary maximum is found at any 
pressure. Fig. 2 (b) are the results for the liquid with a 
critical point and a coexistence line between a LDL phase 
and a HDL phase. Note how the specific heat, for pres- 
sures Pv /e > 0.89, has two maxima. The characteristics 
of the secondary peaks are the following: they are located 
always at a temperature (T*) larger than the nucleation 
temperature, they are non-divergent, and the tempera- 
ture (T*) is pressure independent. These findings agree 
with our hypothesis previously presented and suggest 
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FIG. 2: Constant pressure specific heat vs. tem- 
perature for different pressures Pv /e =0.54 (thiner 
line),0.74,0.84,0.89,0.94,0.955,0.97 (thicker line) and for (a) 
J ff /e = (singularity free scenario) and (b) J a /e — 5 (liquid- 
liquid critical point scenario). Note how we find secondary 
maxima only in (b) and for Pv /e > 0.89, indicating that the 
critical pressure of the model must be close to Pv D /e = 0.89. 
The relative values of the maxima of the peaks are system- 
size dependent expect for the secondary peaks indicated in 
(b). Absolute specific heat values are normalized to Co, the 
value of the secondary maxima found in (b). 



that there is a critical point at a pressure Pv /e ~ 0.89 
for liquid (b). 

There is an important comment about the values of the 
main maxima of the specific heat shown in Fig. 2. Since 
we are dealing with an small system with finite size ef- 
fects the concrete value of the maxima are not represen- 
tative. For example, for the isothermal compressibility, 
it is found that maxima are proportional to the number 
of molecules in the first order phase transition, and that 
they scale as a power of the number of molecules when 
the transition is at the critical point (second order phase 
transition) In any case, they are not measurable in a 
real experiment due to the homogeneous crystallization 
of water. On the contrary, the secondary maxima shown 
on Fig. 2 (b) are clearly non-divergent because they are 
size independent (not shown). 

This non-divergent secondary peak found for P > P c 
could be measured in a real, high pressure, water ex- 
periment, since it is located at temperatures above the 
nucleation temperature. The scheme of the concrete ex- 
periment to be performed is presented on Fig. 3. Calori- 
metric measurements on water should be performed for 
pressures ranging from atmospheric pressure to approx- 
imately P — 2 kbar in the supercooled region. Until 
now, calorimetric measurements in bulk water have been 
mainly performed at low pressures (Icj | and no maxi- 
mum has been found down to T — 240K. However, fol- 




FIG. 3: Schematic representation of the experiment proposed 
to test the existence of a second liquid-liquid critical point at 
P c ~ 1.7 kbar and T c ~ 200K (black circle) . Constant pres- 
sure specific heat measurements should be performed by cool- 
ing down water until reaching the homogeneous nucleation 
temperature (squares). For pressures smaller than ~ 1.7 kbar 
the specific heat should present the usual constant power- 
law increase when cooling in the deeply supercooled region. 
However, for larger pressures (in the range form 1.7kbar to 
2kbar), there should be a previous non-divergent maximum 
in T — T„ and close to the temperature of 200K indepen- 
dently of the pressure applied. This behavior for the specific 
heat has been schematically represented on the figure (thick 
continuous lines). Dotted line is the proposed coexistence 
line. 



lowing our hypothesis, if there is a liquid-liquid critical 
point, calorimetric measurements should present a non- 
divergent maxima for values of the pressure P > P c ~ 
1.7kbar [ll| at a temperature T* above the homogeneous 
nucleation temperature (which is located between 200K 
and 180K depending on the value of the pressure ap- 
plied). The value of T* should be almost pressure inde- 
pendent and close to the hypothesized critical tempera- 
ture T c ~ 2QQK This secondary maximum should 



not show up in other quantities such as the compressibil- 
ity because the intra-molecular term we are considering 
for the water molecule does not imply any change in vol- 
ume. On the other hand, if the real scenario in water is 
singularity-free there should be no maxima in the specific 
heat. 

To conclude, the liquid-liquid critical point and the ex- 
istence of an internal correlation in the water molecule 
could be strongly supported by the existence of non- 
divergent maxima in the water specific heat at high pres- 
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sures. Of course, an important question is: in case these 
experiments are performed and no maxima are found for 
any pressure, does it immediately implies that there is no 
second critical point in water? The answer is no, because 
there is also the possibility of the existence of a liquid- 
liquid critical point in water that is not related to an 
internal correlation inside the water molecule due to the 
existence of an intra-molecular term in the Hamiltonian. 
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